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Abstract. We study theoretically the efficiency of an asymmetric superconducting 
quantum interference device (SQUID) which is constructed as a loop with three 
capacitively and resistively shunted Josephson junctions. Two junctions are placed 
in series in one arm and the remaining one is located in the other arm. The SQUID 
is threaded by an external magnetic flux and driven by an external current of both 
constant (dc) and time periodic (ac) components. This system acts as a nonequilibrium 
ratchet for the dc voltage across the SQUID with the external current as a source of 
energy. We analyze the power delivered by the external current and find that it strongly 
depends on thermal noise and the external magnetic flux. We explore a space of the 
system parameters to reveal a set for which the SQUID efficiency is globally maximal. 
We detect the intriguing feature of the thermal noise enhanced efficiency and show 
how the efficiency of the device can be tuned by tailoring the external magnetic flux. 
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1. Introduction 

The SQUID is the most sensitive instrument which is capable of detecting and measuring 
even extremely small magnetic helds. It has been used successfully not only for 
magnetometry but also for voltage and current measurements. Its applications go far 
beyond the research laboratories often into commercial apparatus exploited in metrology, 
geophysics and medicine, see the reviews HIE]. The SQUID has been the topic of various 
extensive theoretical and experimental studies. Yet, a number of open problems of this 
setup still remain to be resolved. A prominent example may be the efficiency of the 
SQUID as a thermodynamical machine converting the input energy into its other forms. 
It is the subject of this paper. 

We study an asymmetric SQUID driven by an external current and analyze the 
charge transport and voltage induced across the device. The asymmetric SQUID is 
modeled as a ratchet far from equilibrium, i.e. as a classical Brownian particle moving 
in a spatially periodic potential with broken reflection symmetry and driven by a time- 
dependent force. In this mechanical analogy, the voltage across the SQUID corresponds 
to the particle velocity. The most basic measure for characterizing the motion of the 
Brownian particle is its long time average velocity (v). However, alone it does not give 
any information on quality of transport process. Is it effective or ineffective? To answer 
this question, we need to consider its other attributes. One of them are fluctuations 
of the velocity around its average value which in the long time regime are represented 
by the variance Then, typically the instantaneous velocity v(t) takes 

values within the interval of standard deviation, v{t) G [(v) — (v) + ay]. Note that 

if fluctuations are large, i.e. if > |(n)|, then it is possible that the particle moves 
for some time in the direction opposite to its average velocity (n), spread of velocities 
is large and overall transport is not effective. The next feature which is important in 
answering the question about the quality of transport phenomenon is related to the 
ratio of energy input into the system and its energetic output. How much of the energy 
input is converted into directed motion of the particle and how much of it is wasted 
by spreading out into environment and dissipated as heat? A proper quantifier to 
characterize this aspect of transport is the efficiency of the system. 

By using the correspondence between the SQUID and the mechanical ratchet 
system, we study three measures for evaluation of transport quality: average voltage, its 
fluctuations and the efficiency of the SQUID. In the previous paper [3] we analyzed the 
average voltage in this setup for wide parameter regimes: covering the overdamped to 
moderate damping regime up to its fully underdamped regime. We found the intriguing 
features of a negative absolute and differential conductance, repeated voltage reversals, 
noise induced voltage reversals and solely thermal noise-induced ratchet voltage. We 
identihed a set of parameters for which the ratchet effect is most pronounced and showed 
how the direction of transport can be controlled by tailoring the external magnetic 
flux. The main emphasis of that work laid on formulating and exploring conditions 
that are necessary for the generation and control of transport 0 0, its direction. 
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Figure 1. Schematic asymmetric SQUID composed of three Josephson junctions and 
the equivalent circuit composed of two junctions. The Josephson phase difference is 
‘^1 = + the externally applied current is I, the current through the left and right 

arms is R and R, respectively. The external magnetic flux is <I>e and the instantaneous 
voltage across the SQUID is V = V{t). The long time average voltage (U) across the 
SQUID is expressed by the relation (U) = h{ipi) f2e = hlfp-Rl^e. 


magnitude as well as its dependence on system parameters. However, apart from these 
well investigated questions other important features concerning the quality of transport 
iBIH] have remained unanswered. Therefore in this paper we concentrate on this topic 
and connection between the directed transport expressed in terms of the dc voltage, its 
fluctuation characteristics and energetics of the SQUID. 

Theoretical aspects considered in the paper concern not only our specific SQUID 
ratchet but a much wider class of systems and problems. There are many experiments 
on a number of ratchet systems |9], in particular superconducting ratchets [10], a part 
of which can be controlled by an external magnetic held as well as 

theoretical studies of such systems driven by harmonic and biharmonic external currents 
[IS] dSl El UHl DSj. However, the efficiency of transport has not been analyzed in the 
above-cited papers. 

The structure of the paper is as follows. In Sec. H, we recall the model of a 
SQUID rocking ratchet which is composed of three resistively and capacitively shunted 
Josephson junctions. In Sec. HI we dehne mean values of arbitrary state functions in 
the long time regime. Then in Secs. IV-VI, the quantities characterizing the quality 
of the transport such as the voltage huctuations, the energy balance and the (Stokes) 
efficiency are introduced, respectively. In Sec. [^we elaborate on key aspects of transport 
efficiency in the system: starting from the power delivered by the externally applied 
current, covering the tailoring of the Stokes efficiency of the device, up to presentation 
of the regime for which thermal noise enhances the efficiency and discussion about the 
impact of variation of the external magnetic flux on the efficiency of the SQUID. Finally, 
the last section provides a summary. 
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The asymmetric SQUID [TUI ED ED ED El] is presented in Fig. 1. It is a loop with two 
resistively and capacitively shunted Josephson junctions [25] in the left arm and one in 
the right arm. The crosses denote the junctions and ipk = {k = n, d, 1, 2) are the 

phase differences across them. Each junction is characterized by the capacitance (7^, 
resistance and critical Josephson supercurrent respectively. 

The SQUID is threaded by an external magnetic ffux $e and driven by an external 
current I = I{t) which is composed of the static dc current Jq and the ac component of 
amplitude A and angular frequency D, namely 

I{t) = Iq + Acos{Ilt). (1) 

To reduce a number of parameters of the model, we consider a special case when two 
junctions in the left arm are identical, i.e. Ju = Jd = Ji,Ru = Rd = Ri/‘^,Cu = Cd = 
2Ci. In some regimes [3], two junctions in a series can be considered as one for which 
the supercurrent-phase relation takes the form Jisin ((pi/2), where ipi = (pu + (pd- This 
result is also derived in Ref. [26| for an effective double-well structure described in terms 
of a double-barrier potential (cf. Eq. (23) therein). 

The total magnetic ffux 4) piercing the loop is a sum of the external ffux and 
the ffux due to the ffow of currents. 


$ = <he Li, 


( 2 ) 


where L is the loop inductance and i = i{t) is the circulating current which tends to 
screen the magnetic ffux. In the ’’dispersive” operating mode of the SQUID 1271 , i.e. 
when the condition \Li\ -C $0 holds true (<Fo = h/2e is the ffux quantum), the phase 
(f = (fi obeys the Stewart-McCumber type equation of the form |3| 


where the effective supercurrent J{p) reads 



( 3 ) 


J{(p) = + J2sm{(p + ^e)- (4) 

The parameters are: C = Ci + C 2 , R~^ = Rf^ + Rf^, is the Boltzmann constant, T 
is temperature of the system and <Fe = 27r$e/‘ho is the dimensionless external magnetic 
ffux. Thermal ffuctuations are modeled by J-correlated Gaussian white noise f{t) of 
zero mean and unit intensity 


m) = o, = (s) 

The Stewart-McCumber equation (|^ has the form of a Langevin equation and describes 
a non-Markovian stochastic process for the phase p. In the extended space {p,(p}, it 
has the Markovian property and all well known methods can be applied to analyze it. 

Eq. ([^ can be interpreted in the framework of a model of a classical Brownian 
particle. It helps develop the intuition and interpretation. In the one-to-one 
correspondence, the particle position x translates to the phase p, the particle velocity 
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V = X to the voltage V cc (p, the conservative force to the supercurrent J{(p), the 
external force to the current I{t), the mass m to the capacitance m oc C and the 
friction coefficient 7 to the normal conductance 7 oc G = 1/i?. It is important to 
note that the friction 7 is not proportional to the normal resistance R (as one could 
expect in the case of electrical circuits) but to the inverse of R. The reason is that 
plasma oscillations of the junction are more damped if more normal electrons couple 
to the oscillating condensate (i.e. when G is greater). The voltage accelerates normal 
electrons and their kinetic energy is dissipated into heat. Thus the plasma oscillations 
converts into heat with a rate proportional to the conductance G [28] . 


3. Asymptotic mean values 


The main characteristics of the system are the current-voltage curves in the long time 
regime. It can be shown that for the external current ([^, they can be extracted from 
the relation for the averaged voltage (y) developed across the SQUID [3] 




( 6 ) 


where pi can be obtained from the Stewart-McCumber equation ([^ for the phase 
p = pi. The above relation holds true when the average is only over the period 
T = 2tt/II of the ac current. It is also valid in the long time regime when the averaging 
is performed over initial conditions and all realizations of thermal fluctuations. Below 
we present a precise definition of this procedure. 

It is convenient to rewrite Eq. ([^ for the phase p and the voltage V = V{t) = 
{h/2e)p in the Ito form 


2 e 

dp = —Vdt, 
n 






dt+t 


2kBT 


R 


dW{t), 


( 7 ) 

( 8 ) 


where dW{t) = f{t)dt is the differential of the Wiener process of zero mean and the 
second moment {dW{t)dW{t)) = dt. The pair {p, V} form a Markovian process and its 
probability density P = P{p, V,t) obeys the Fokker-Planck equation 


dP 2e,, d ^ Id 

— =- V—P H- 

dt h dp G dV 


+ J{p) 


lit) 


P + 


knT d^ 


-P 


( 9 ) 


RG'^ dV^ 

with the initial condition Pip, V, 0) = pip, V), where a given probability density pip, V) 
describes the initial distribution of the phase (p( 0 ) and voltage U( 0 ). 

For any state function fip,V), its mean value {fip,V))t at time t is calculated 
from the relation 

p27T poo 

{fip,V))t= dp dV fip,V)Pip,V,f). (10) 

-70 J —00 

Because the system is driven by the time-periodic current lit), the probability density 
Pip,V,t) approaches for long time the asymptotic periodic form Pasip,V,f), namely 
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Pas{y:>,V,t) = ( 11 ) 

n=—oo 

where the Fourier coefficients Wn{y^,V) are solutions of the differential equations 
obtained from the Fokker-Planck Eq. §. The time-dependent asymptotic mean value 

^27r f*oo 

{fi^:V))r= dV f{^,V)Pas{^,VC) ( 12 ) 

Jo J —CO 

is also a periodic function of time. If we are interested in its time-independent form, 
the time averaging over the additional period T = 2tt/II of the ac current has to be 
performed 

Q ^t+27r/r2 

{f{^,v)) = - {f{ip,v)r:du 

Q pt-\- 27 v/Q 

= / {f{^,V))udu. (13) 

In a particular case, when f{(p,V) = V, we get the time-independent asymptotic mean 
voltage (V). Similarly, when f{(p,V) = V^{k = 2,3,...), we obtain the stationary 
statistical moments of the voltage {V^). 


4. Fluctuations of voltage 


The asymptotic average voltage (V) calculated according to the prescription (13) is 
the most important transport characteristics of the system. The magnitude of the 
instantaneous voltage V{t) can be much larger than its mean value. Moreover, the 
fluctuations of voltage in the long time regime can also be large. They are described by 
the voltage variance 


4 = {v^) - {v)\ 

The voltage typically ranges within the interval of several standard deviations 


(14) 


V{f) e[{V)-nav-,{V)+nav\, n = l,2,... (15) 

If the standard deviation ay is large, i.e. when ay > |(1^)|, the voltage V{t) can spread 
far from its average value and even assume the sign opposite to it. It is a case for protein 
motors in biological cells where the instantaneous velocity changes direction very rapidly 
and its absolute value is several orders of magnitude larger than the average velocity 

ra. 


5. Energetics of the SQUID 

The SQUID is a device which converts input energy into its other forms. It is provided 
by the external current I(t) and the energy flow is determined by the equation of motion 



Efficiency of the SQUID Ratchet Driven by External Current 


7 


(|^. In the mechanical interpretation, the kinetic energy of the particle corresponds to 
the energy stored in the system of capacitance C, namely 


Ec = Ec{V{t)) = -CV^it). 


(16) 


The particle potential energy translates to the Josephson energy accumulated in the 
junction when the supercurrent flows through it 

Ej = Ej{(p{t)) = [ J{ip{u))V{u)du 


^ ^Jo ’^(7>iu))(p{u)du= — 


h , ^ 

= ~ ^ ^ 2 cos 


ip{t) 


J{(j))d(j) 

+ J 2 cos [93(f) + l>e] 


const. 


(17) 


The sum 


E = Ec + Ej (18) 

is the total energy of the system. Its balance can be obtained from Eqs. (0-(§. For 
this purpose we apply the Ito differential calculus to both functions Ec{V) and Ej{ip) 


dEc — 






J{^)V + I{t)V + 


keT 

RC 


dt T 


2,k]^T 


R 


dEj = 


dE 


d^p 


dtp = J{ip)Vdt. 


VdW{t), (19) 

( 20 ) 


Next, for both sides of the above equations, we calculate the mean values. Exploiting 
the Ito martingale property we hnd the average value of the term {VdW{t))t = 0 and 
obtain the energy balance equation in the form [HI E2] 




RC 


1 


= -^[{V^)t-{V^)e,]+I{t){V)c 


( 21 ) 

( 22 ) 


where (■)* denotes a mean value at time t according to the prescription (10). In the 


right hand side of this equation there are three components, each of them is related to 
the separate process responsible for the energy change. Let us point out that the hrst 


term in Eq. (21) is always negative whereas the third is positive. The former describes 


the rate of energy loss due to dissipation and the latter refers to the energy provided 
by thermal equilibrium fluctuations. According to the equipartition theorem, in the 
thermodynamical equilibrium, when I{t) = 0, the relation {CV‘^l2)eq = ksT/2 holds 


true. It is utilized in Eq. (21) to get 


The second term in (21) characterizes the 


change of energy caused by the external current I{f). 

Now, we perform the hnal averaging of Eq. ([2l|) according to the prescription (13) 


Since in the long time regime the average values are periodic function of time, the left 
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hand side vanishes 

/ t+2'tT/Q 1 ^ 

-{E)^du = -[{V\t + 2n/Il))-{V\t))] 

+ [{Ej{^{t + 27r/(]))) - {Ej{^{t)))] 

= 0. (23) 

As a conseqnence, in the stationary regime the mean power Vin delivered to the system 
by the external cnrrent I{t) over the period T is expressed by the relation [32] 

= (Ht)V) = ^ [{V^) - {V\,] . (24) 

From the above eqnation it follows that the amonnt of the energy inpnt to the SQUID 
from the external driving I{t) depends not only on the cnrrent itself (i.e. on /o,A,D) 
bnt also on properties and parameters of the device: its temperatnre T, the resistance 
R and the capacitance C . In contrast, the energy snpplied by thermal flnctnations does 
not depend on the external current but only on T, R and C. 


6. Efficiency of the SQUID 


A generic definition of the efficiency of a device converting energy is a ratio between the 
output (work, power) and the input (power) energy. Depending on the choice of input 
and output, different definitions of the efficiency characterize various aspects of energy 
conversion in the device. To explain the problem, we use the mechanical interpretation 
of (j^. Then the average voltage (U) corresponds to the velocity of the Brownian particle 
or the Brownian motor. The thermodynamic efficiency is defined as a ratio of the work 
done by the motor to the energy input. If the particle is working against a constant 
force (load) i^ then in stationary state the efficiency is defined as follows 

io{V) 


Vo = 




(26) 


In the considered case, there is no a load and the Brownian motor does not transport 
external objects. It works against the friction ’’force” 


Ir =- 

" 2eR^ 


(26) 


(By the way, it has not the unit of Newton but if the above formula is again multiplied 
by the factor h/2e then it has the correct physical unit.) When the external force I{t) 
is switched off, the velocity of the motor is damped to zero and the system tends to 
thermodynamical equilibrium. Because the motor works against the friction force, we 
can utilize its mean value to get another definition of efficiency, namely 

^ {Fn){V) {VY ^ {V? , . 

V,n RV,n {VY-{VYe, {V^ + a^y - ksT/C' ^ ’ 

This quantity is called the Stokes efficiency [SS] [3ll |3l] . Let us note that it depends 
explicitly on mass C of the Brownian particle and and only implicitly on the friction 
coefficient R via the Langevin equation ([^. 
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It should be mentioned that (27) is not the rate of the work done by the motor on 
its surroundings (viscous medium). Moreover, it is not a mean power Vr to overcome 
the friction force which correct form reads 

(V^) 


Vr = {FrV) = 


R 


(28) 


However, this expression cannot be put as a numerator in the definition of the efficiency 
because there are regimes where the mean velocity is extremely small (numerically zero), 
(y) Tn 0 but y^) 7 ^ 0 and the efficiency could be large even though the particle does 
not move on average in one direction. It is the main reason why the Stokes efficiency is 
more adequate is such cases as considered in the paper. 

Another possible definition of the efficiency is based on the remark that what we 
observe in the long time regime is the average velocity. Therefore we can introduce 
’’kinetic power” as the ’’kinetic energy” of the particle per the period T, 

cyf 


Vk = 


2T 


(29) 


One should note that it is not exactly proper definition of the kinetic power as it should 
be proportional to (V^) instead of (H)^. However, we replaced it with the latter for the 
reason explained before. Nevertheless, it is still a measure of performance of the motor. 
If the average velocity increases then Vk also grows. We can insert it as a numerator in 
(25) and then we get the kinetic efficiency 
Vk RC (H)2 


rjk = 


Vi, 


2T (H2) - 


2T 


(30) 


This quantifier can be used only when the time-periodic force is switched on. Then 
qualitatively, it is similar to the Stokes efficiency. However, the dependence on the mass 
C, the friction coefficient R and the period T is different. Both the Stokes efficiency 
and the kinetic efficiency are consistent with our intuition: a decrease of fluctuations ay 
leads to a smaller input power and hence to an increase of the efficiency. Consequently, 
the transport is optimized in regimes that maximize the directed velocity and minimize 
its fluctuations. Because the kinetic efficiency is proportional to the Stokes efficiency, 
below we analyze only the last one. 


7. The results 

7.1. Dimensionless model 

There are several dimensionless forms of Eq. (|^ in dependence of the choice of 
scaling time. In this system there are four characteristic frequencies: plasma frequency 
oOp = 2eJi/hC, the characteristic frequency of the junction ooc = 2eRJi/h, the frequency 
Ur = 1/RC related to the relaxation time and the frequency 12 of the ac current. There 
are three independent characteristic time scales related to these frequencies (note that 
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j = 0 - 

=7r/2 -- 
—7r/2 —• 



Figure 2. The potential (33) for the symmetric case j = 0 (solid red line) in 
comparison with the ratchet potential for j = 1/2 and two values of the external 
magnetic flux l>e = 7r/2 (dashed green line) and $e = —r 12 (dotted blue line). 


ujp = UcOJr)- Here, we follow pni and define the new phase x and the dimensionless time 
t as 


X = 


</? + vr 


S = —, Tc = 


h 


2 ' ’ Tc " eRJi 

Then ([^ takes the dimensionless form 

Cx{s) + i;(s) = —U'{x{s)) + F + acos{ujs) + \/2D f{s), 


(31) 

(32) 


where the dot and prime denotes a differentiation over the dimensionless time s and the 
phase X, respectively. We introdnced a spatially periodic potential U{x) of period 27r of 
the following form 121 

U{x) = — sin(x) — ^ sin(2x + l>e — 7r/2). (33) 

This potential is reflection-symmetric if there exists xq snch that U (xq + x) = U (xq — x) 
for any x. If j ^ 0, it is generally asymmetric and its reflection symmetry is broken, 
see Fig. We classify this characteristics as a ratchet potential. However, even for 
j 7 ^ 0 there are certain valnes of the external flnx <Fe for which it is still symmetric. 
The dimensionless capacitance C is the ratio between two characteristic time scales 
C = Tr/Tc^ where the relaxation time is Tr = RC. Other re-scaled parameters are 
j = J 2 /J 1 , F = Iq/Ji, a = A/Ji and oj = Qtc. The rescaled zero-mean Ganssian 
white noise ^(s) has the anto-correlation fnnction {f{s)f{u)) = 6{s — u) and its intensity 
D = eksT/hJi is the qnotient of the thermal and the Josephson conpling energy. The 
dimensionless voltage v{t) = x(s) = Vif)/RJi and therefore the physical average voltage 
(y) is given by the relation 

{V) = RR{v). (34) 


In particular, after such a scaling procedure the dimensionless input power Pin is 
expressed as 


Pin = {t>") - Die 


(35) 
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and consequently, the Stokes efficiency reads 

(n)2 + a2 - D/C {v^)-D/C ^ ^ 

The key feature for the occurrence of the directed transport {v) 7 ^ 0 is the symmetry 
breaking. This is the case when either the dc current F 7 ^ 0 or the reflection symmetry 
of the potential U{x) is broken. 

The system described by (32) becomes deterministic when the thermal noise 
intensity D is set to zero. Even in this case it exhibits complex dynamics including 
chaotic regimes [351136] . The application of noise generally smooths out its characteristic 


response function. There are two classes of states of the driven system dynamics: the 
locked states, in which the phase stays inside hnite number of potential wells and the 
running states for which it runs over the potential barriers. The latter are crucial for 
the occurrence of the transport. They can be either chaotic (diffusive) or regular. 


1.2. Details of simulations 

The Fokker-Planck equation ([^ corresponding to the Langevin equation ([^ cannot 
be solved by use of closed analytical forms. Therefore, in order to obtain the relevant 
transport characteristics we have to resort to comprehensive numerical simulations of the 
driven Langevin dynamics. We have integrated (32) by employing a weak version of the 
stochastic second order predictor corrector algorithm 1371 with a time step typically set 
to about 10“^-27r/a;. Since (32) is a second-order differential equation, we have to specify 
two initial conditions a;(0) and i;(0). Moreover, because for some regimes the system may 
be non ergodic in order to avoid the dependence of the presented results on the specihc 
selection of initial conditions we have chosen phases x( 0 ) and dimensionless voltages 
i;(0) equally distributed over interval [ 0 , 27r] and [—2,2], respectively. All quantities of 
interest were ensemble-averaged over 10 ^ — 10 ^ different trajectories which evolved over 
10^ — 10^ periods of the external ac driving. Numerical calculations were done by use 
of a CUBA environment implemented on a modern desktop GPU. This scheme allowed 
for a speed-up of a factor of the order 10 ^ times as compared to a common present-day 
CPU method [38l |39]. Part of our so obtained results are presented next. 


7.3. Power delivered by external current 

Let us begin analysis of the SQUID efficiency by looking at the power Pin delivered by 
the external current I{t). Notably, it depends implicitly not only on the parameters 
of the applied external current {F, a, u) but also on the quantities characterizing the 
device like the capacitance C. We have found that generally the input power ( [^ 
tends to increase for larger values of the dc current F and ac driving amplitude a. The 
dependence on the frequency u is more complex. However, in most cases Pin is relatively 
large when u is small. It is because very fast oscillation of the driving current cannot 
induce neither the average voltage (n) nor (n^). In panel (a) of Fig. we show the 
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Figure 3. The power (241 delivered by the external current I(t) is presented in the 
dimensionless form (35) as a function of the dimensionless capacitance C, the thermal 
noise intensity D and the external magnetic flux <i>e in panels (a)-(c), respectively. 
Other parameters read: a = 1.9, uj = 0.6, j = 1/2. In panel (a): D = 9.7 ■ 10“® and 
$e = r-/2. In panel (b): C = 0.645 and $e = r'/2. In panel (c): C = 0.645 and 
£> = 9.7-10-5. 


representative dependence of the input power on the dimensionless capacitance C 
of the SQUID. One can observe that Pin is maximal for the overdamped or close to 
damped regime and decreases when C grows. Since in the mechanical framework the 
capacitance C translates to the mass of the Brownian particle it is intuitively clear that 
when the inertial term becomes large then the device needs more power to response 
equivalently. Perhaps the most surprising is the fact that Pin depends explicitly on the 
thermal noise intensity D, i.e. on temperature of the system. Typically, it decreases 
for increasing D. However, there are also regimes for which Pin is enhanced by thermal 
noise. In panel (b) of Fig. [^we exemplify this situation. Indeed, for a wide interval of 
temperature the input power is almost monotonically increasing function of the noise 
intensity D. Finally, the influence of the constant external magnetic flux on Pin is 
depicted in the last panel. It is remarkable that one can tune the input power Pin by 
changing the external magnetic flux. The reader should note that for the presented 
regime it is maximal when = 0, i.e. potential U{x) is reflection symmetric. In such a 
case there is no average voltage drop (v) = 0 across the device when additionally the dc 
current F vanishes. It follows that large input power Pin does not necessarily translate 
into the efficient directed transport. 


l.f. Tailoring Stokes efficiency 


The system described by Eq. (32) has a 7-dimensional parameter space 
{C,a,u;,F,j,^g,D}. We set F = 0 and check how it depends on the remaining 
system parameters. We limit our considerations to positive a because the system 
(32) is symmetric under changes of sign of a. Depending on the magnitude of the 


dimensionless capacitance C of the device it can operate in three distinct regimes: 
overdamped (C —)• 0), damped (moderate C) and underdamped (C —)• oo). We 
note that the conditions that are necessary for the generation and control of the 
direction of transport have been extensively studied in these regimes in our recent 
work [3]. Since very fast oscillation of the driving current cannot induce the average 
voltage (v) it is sufficient to limit our considerations to low and moderate ac driving 
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Figure 4. The Stokes efficiency ys defined by Eq. (27) in the parameter plane {a,a;} 
of the ac current for three distinct regimes: overdamped (C = 0.2), moderate damping 
{C = 2) and underdamped {C = 10) in panels (a), (b) and (c), respectively. The 
remaining parameters are: F = 0, D = 10“®, j = 1/2 and $e = r/2. 


frequencies u. We have performed scans of the following area of the parameter space 
C X a X u E [0.1; 10] x [0; 10] x [0.1; 1] at a resolution of 200 points per dimension to 
determine the general behavior of the system. The results are depicted in Fig. 

We can see that regardless of the regime in which the device operates its Stokes 
efficiency rjs is zero or negligibly small for a < 1 and high frequencies u. This is due 
to the fact that the rocking mechanism is either too weak or too fast to induce hnite 
average voltage (v). The areas of non-zero efficiency rjs have a striped structure. For a 
given amplitude a, the ratchet behavior generally tends to disappear as the frequency 
uj grows. On the other hand, for a given frequency, there is optimum amplitude a that 
maximize the Stokes efficiency. The increase of the capacitance C causes blurring of the 
regions for which the efficiency is nonzero. Moreover, this tendency is often accompanied 
by its reduction. Consequently, the studied device operates best in the overdamped or 
close to damped regimes. 


7.5. Optimal regime 


We have explored the parameter space of the system (32) and we have been able to detect 
a regime for which the efficiency ps is globally maximal. It is in the vicinity of the point 
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Figure 5. Optimal regime for the Stokes efficiency of the transport occurring in (321. 
The dependence of the average voltage {v), its fluctuations and finally the Stokes 
efficiency ys on the static dc-bias F, frequency uj, amplitude a, capacitance C and 
thermal noise intensity D is presented in panels (a)-(o). Parameters are: F = 0, 
UJ = 0.406, a = 1.55, C = 0.496, D = 10-^ j = 1/2 and = tt/2. 


{C,a, o;} = {0.496,1.55,0.406}. All relevant transport characteristics, i.e. the average 
voltage (v), its flnctnations and the efficiency rjs corresponding to neighborhood of 
this set of parameters are presented in Fig. 5. 

Let ns begin with the dependence of the transport characteristics on the dc cnrrent 
F. It is depicted in Fig. |^a)-(c) for small valnes of F G (—0.5, 0.5). Panel (a) presents 
the cnrrent-voltage cnrve. In the low temperatnre limit {D = 10“^) the average voltage 
is almost qnantized at valnes nu, n = 0, ±1,... For a symmetric potential, these plateans 
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correspond to standard Shapiro steps m However, in our case also steps at half 
integer multiplies of u can be observed. This is due to the deviation of U{x) from a 
simple sin a: form, which is the sole case for which steps lie only at integer values of 
nuj [20]. However, in both the symmetric and asymmetric cases a proper amount of 
noise is sufficient to wipe out their evident structure |in|. This Shapiro-like current- 
voltage curve is characteristic for the device operating in the low temperature limit of 
overdamped or damped regimes. Panel (b) of the same Fig. 5 presents the dependence of 
the voltage fluctuations on the dc current F. It is rather complicated non-linear and 
non-monotonic function of F without any immediately obvious relation to the average 
voltage of panel (a). However, the most important observation is that the voltage 
fluctuations are minimal for F = 0. This fact is of fundamental importance for the 
influence of F on rjs- In Fig. |^c) we can see that rjs is locally maximal for F = 0. For 
large values of F (not shown here) the mean voltage is an almost linear function of F 
and the efficiency approaches the value 1. It support the statement in Ref. [31] that 
when the Stokes efficiency is close to 1, the driving resembles the constant force. 

The role of the frequency u of the ac current is illustrated in Fig. [^d)-(f). Panel 
(d) presents the dependence of the average voltage (v) on u. In the adiabatic limit 
a; —0 it undergoes rapid oscillations. This behavior has its reflection in the influence 
of the frequency on the voltage fluctuations (see Fig. |^e)). According to the 
previous statement very fast oscillations of the driving current cannot induce the non¬ 
zero average voltage. Therefore for sufficiently high frequency there is no transport and 
as a consequence the Stokes efficiency rjs is zero. However, a strong peak of efficiency 
is observed for moderate value of ca = 0.406. It is associated with the fact that for this 
frequency the average voltage (v) is maximal and simultaneously its fluctuations are 
minimal. 

The impact of the amplitude a of the ac current is shown in Fig. [^g)-(i), 
respectively. In particular, the resonance-like behavior is observed in the dependence 
of the average voltage (v) on the amplitude a (see Fig. [^g)). Apart from two clearly 
visible peaks there is almost imperceptible small directed transport. This fact has 
critical impact on the functional dependence of the efficiency rjs. It is proportional to 
(n)^ so it vanishes too when the device response is zero. The influence of variation of the 
amplitude a on the voltage fluctuations is depicted in panel (h). It is almost linearly 
increasing function of a. Only one evident deviation from this trend can be noted, i.e. 
a local minimum around a = 1.55 corresponding to the hrst high peak in Fig. [^g). It 
should be stressed that there is no any contradiction with the dependence of the average 
voltage since (v) = 0 does not necessarily mean (n^) = 0 and therefore can at the 
same time assume nonzero value. 

The dependence of all relevant transport characteristics on the capacitance C is 
very complicated as shown in Fig. [^j)-(l). The hrst panel of this group shows the 
average voltage versus the capacitance C. We note the important feature of the voltage 
reversal [2S1 SI] : starting from zero, the voltage changes its sign from positive to negative 
and again in the opposite direction as C grows. Therefore the capacitance can serve 
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Figure 6. Noise enhanced Stokes efficiency. The dependence of the average voltage 
(n), its fluctuations cr„ and finally the efficiency ys on thermal noise intensity D is 
presented in panels (a)-(c), respectively. Parameters read: a = 1.899, uj = 0.403, 
C = 6, j = 1/2 and = 7r/2. 


as a parameter to manipulate the direction of transport processes. The efficiency rjs is 
maximal close to the border between the overdamped and damped regimes. It is (almost) 
zero in the underdamped limit which corresponds to large capacitance C* —?• oo. This 
is a consequence of the fact that for this regime the average voltage (v) vanishes or is 
negligibly small. 

The last three panels of Fig. [^depict the influence of the thermal noise intensity D 
on all previously studied quantities. An increase of the noise intensity D leads to both 
monotonic decrease of the induced average voltage (v) and increase of its fluctuations 

Consequently, the efficiency is the best in the low temperature regime when the 
deterministic dynamics of the system (32) plays a crucial role. 


7.6. Noise enhanced Stokes efficiency 

We have found the opposite scenario when thermal noise enhances the efficiency. This 
perhaps surprising effect is exemplified in Fig. Panel (c) presents the dependence of 
the efficiency rjs on D. Evidently, in some intervals of D, the increase of temperature 
causes the increase of t]s. There is also optimal value of temperature or equivalently 
the thermal noise intensity D ^ 0.0004 for which the efficiency takes its maximum. 
Moreover, in this case the ratchet mechanism is solely activated by thermal equilibrium 
fluctuations as for low noise intensity no rectification can be observed. This statement is 
con&med in the functional dependence of the average voltage (v) which is presented in 
panel (a). It is also remarkable that in this regime an increase of thermal noise intensity 
D leads to a decrease of voltage fluctuations cx^. 


7.7. Impact of external magnetic flux 

As it was shown before, the efficiency can be tuned in several ways. However, it seems 
that from the experimental point of view the simplest method is to vary the external 
constant magnetic flux <Fe. The dependence of the average voltage (n), its fluctuations 
ay and the efficiency on the external magnetic flux $e in the previously presented regime 
for which thermal noise induces the ratchet effect (cf. Fig. is depicted in Fig. From 
the symmetry considerations of (32) it follows that for an arbitrary integer number n, the 
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Figure 7 . Impact of the external magnetic flux on the average voltage (v), its 
fluctuations cr„ and Anally the efflciency ys is presented in panels (a)-(c) for the 
noise intensity D = 0.0004 which corresponds to the maximum in Fig. 6(c). Other 
parameters are the same as in Fig. 



Figure 8. The Stokes eSiciency of the rocked SQUID in the parameter plane {$eO}- 
The left panel corresponds to the optimal regime as in Fig. and the remaining 
parameters are the same as there. The right panel presents the regime depicted in Fig. 
I^with D — 0.0004. 


transformation —)■ 27rn —<l>e reverses the sign of the average voltage (n) —>■ —(n). This 
fact can be directly observed in panel (a) of Fig. However, this is not the case for the 
voltage flnctnations where they are symmetric aronnd $e = 0. A carefnl inspection 
of panel (b) reveals that one can reduce magnitude of by nearly two times just by 
correct adjustment of the external magnetic flux. This fact has further consequences in 
the dependence of the efficiency which is depicted in panel (c). It can be slightly tuned 
by a small variation of the external magnetic flux. 

In Fig. 1^ we present how the Stokes efficiency behaves in the parameter plane 
that specihes the form of the spatially periodic potential U{x). For both 
sufficiently small and large j it vanishes completely. One can observe that for a 
given external magnetic flux <he the Stokes efficiency generally tends to increase as 
the parameter j grows. On the contrary, for a given j there is an optimal value of the 
external magnetic flux $e for which the Stokes efficiency is maximal. We note that for 
two presented regimes, the set structure of the non-zero efficiency in the parameters 
plane {<I>e,j} is radically different. The left panel looks like butterfly wings and the 
right is similar to a rocking horse. 
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8. Summary 

In this paper, we comprehensively studied the Stokes efficiency of the asymmetric 
SQUID in the case of the non-zero capacitance of all Josephson junctions and in presence 
of thermal noise. It allowed to analyze transport properties in the system for the 
entire scale of regimes: starting from overdamped, by damped and finally underdamped 
one. We focused on the connection between the directed transport characterized by 
the voltage across the SQUID and its efficiency. In particular, we examined voltage 
fluctuations and energetic performance of the device. We derived the expression for the 
power delivered by externally applied current and discuss its dependence on the system 
parameters. Apart from the expected influence of the current parameters Jq, A and D 
it also depends on the thermal noise intensity D, i.e. on temperature of the system. 

We have found that regions of low efficiency of the SQUID dominates in the 
parameter space. However, we have identified remarkable and distinct regimes of high 
efficiency rjs ~ 0.65. It turns out that the device operates best in the overdamped 
or close to damped regimes. Moreover, with the help of the computational power of 
modern GPU supercomputers we have identified the tailored set of parameters for 
which the efficiency rjs is globally maximal and for this regime we discussed impact 
of variation of almost all system parameters on the relevant transport quantities. In 
particular, it follows that thermal fluctuations often have destructive influence on the 
energetic performance of the device. Moreover, we were able to detect also the regime 
for which thermal noise enhances the efficiency by inducing the large average voltage 
and minimizing its variance. Last but not least, we discussed in detail the impact of 
the external magnetic flux on the performance and effectiveness of the SQUID. 

Our results can readily be experimentally verified with an accessible setup 
consisting of three resistively and capacitively shunted Josephson junctions formed in an 
asymmetric SQUID device. Some partial transport characteristics like voltage have been 
experimentally studied in the overdamped regime [iniEi]. However, the underdamped 
regime has not been tested and the efficiency has not been measured which makes our 
study a challenge for experimentalists. 
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